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Abstract 

Magnetic susceptibility, entropy and specific heat are calculated at the equilib- 
rium points of phase transition to a phase of coexistence of ferromagnetic order 
and superconductivity in a new class of spin-triplet ferromagnetic superconductors. 
The results are discussed in view of application to metallic ferromagnets as UGe2, 
ZrZn 2 , URhGe. 



1. Introduction 

Recently, a spin-triplet superconducting phase has been discovered in a class of new 
ferromagnetic superconductors (UGe2 [1], ZrZn 2 [2], URhGe [3]). In these metallic com- 
pounds the superconductivity exists only in the domain of stability of the ferromagnetic 
order. Moreover, the ferromagnetic order enhances the superconductivity because of the 
unconventional Cooper pairing with spin S =1 [4-7]. A phenomenological model [7] of 
Ginzburg-Landau (GL) type has been used to describe the possible ordered phases and 
the phase diagram of these metallic compounds [8,9]. 

In this paper we present the results of our calculation of magnetic susceptibility, entropy 
and specific heat near the phase transitions in ferromagnetic superconductors with a 
spin-triplet superconductivity. For our aims we shall use the GL free energy introduced 
in Ref. [7] and results published in Refs. [8,9]. 

2. Model and phase diagram 

We consider the GL free energy [7-9] F = J d 3 xf(ip,A4), where 

/ = ^-(D^Y(D^) + a s H 2 + ^| 4 + a f M 2 + ^M 4 + i la M.{^ x V>*) . (1) 
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In Eq. (JTJ), Dj = (V — 2ieAj/hc), and A,- (j = 1,2,3) are the components of the vector 
potential A related with the magnetic induction B = V x A. The complex vector if> = 
(ipx, ip2, "03 ) = {^Pj} is the superconducting order parameter, corresponding to the spin- 
triplet Cooper pairing and M. = {M.j} is the magnetization. The coupling constant 
7o = 47rJ > is given by the ferromagnetic exchange parameter (J > 0). Coefficients 
a s = a S {T — T s ) and a/ = a/(T — Tf) are expressed by the positive parameters a s and 
a/ as well as by the superconducting (T s ) and ferromagnetic (Tf) critical temperatures in 
the decoupled case, when Aiipiipj-mter action is ignored; b > and (3 > as usual. 

We assume that the magnetization Ai is uniform, which is a reliable assumption outside a 
quite close vicinity of the magnetic phase transition but keep the spatial (x—) dependence 
of if). The reason is that the relevant dependence of if) on x is generated by the diamagnetic 
effects arising from the presence of AA and the external magnetic field H [8,9] rather than 
from fluctuations of if) (this effect is extremely small and can be ignored). First term 
in (jT} will be still present even for H = because of the diamagnetic effect created by the 
magnetization Ai = B/An > 0. As we shall investigate the conditions for the occurrence 
of the Meissner phase where if) is uniform, the spatial dependence of if) and, hence, the 
first term in r.h.s. of (0) will be neglected. This approximation should be valid when the 
lower critical field H C \(T) is greater than the equilibrium value of the magnetization Ai 
in the phase of coexistence of superconductivity and ferromagnetism. 

We take an advantage of the symmetry of model and avoid the consideration of equiv- 
alent thermodynamic states that occur as a result of the respective continuous symmetry 
breaking at the phase transition point but have no effect on thermodynamics of the system. 
That is why we shall assume for concreteness of our calculation that the magnetization 
vector is along the z-axis: Ai = (0,0, Ai), where M. > 0. This concrete choice of the 
direction of the magnetization vector does not restrict the generality of the present anal- 
ysis and leads to the same structure of the ordered phases as previously predicted and 
discussed on the basis of general symmetry group considerations in Ref. [7]. 

We find convenient to use the following notations: (fj = ft^tpj, <pj = <fijexp(9j), M = 
p/*M, 7 = 7o/(W /4 , t = a,/Vb, t = a f /y/P, 2 = (0? + 01 + 01), and 6 = (6 2 - 6,). 
We consider the stable homogeneous phases in a zero external magnetic field (H = 0) 
that are described by uniform order parameters Ai and if>. These phases are three: 
(i) the normal phase (0 = M = 0), hereafter referred as "N", (ii) the ferromagnetic 
phase (0 = 0, M > 0), hereafter referred as "FM", and (iii) the phase of coexistence of 
superconductivity and ferromagnetism (0 > 0, M > 0), which will be called "FS"-phase. 
The standard N-FM phase transition occurs at Tf. The N-FM, N-FS and FM-FS phase 
transitions are shown in Fig. 1 [8,9]. The r-axis above the point C coincides with the 
line of the N-FM phase transition of second order. The dashed lines correspond to phase 
transitions of second order: FM-FS transition for t < (— 7 2 /4) and a N-FS transition for 
t > 7 2 /2. The solid lines BC and CA describe first order FM-FS and N-FS transitions, 
respectively. The line of circles BCA defines the borderline of stability and existence of 



2 



r 

2- 

O 

o 

o 

o 

o 

o 

o 

FM 

° 

o 

o 

o 

1- 

'•p 
FS 


N 

,C 

\o 
\. o 
\. o 

\. o 

— , — , — , — , — r-^T-T— -— - r — 7 — -— - 


-1 -0.5 


0.5 1 t 



Figure 1: The phase diagram in the plane (t, r) with two tricritical points (A and B) and 
a triple point C; 7 = 1.2. Second order transition lines are dashed, first order transition 
line is solid. 
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FS. In the domain confined between the line of circles and the dashed line on the left 
of the point B the stability condition for FS is satisfied but the existence condition is 
broken. The phase diagram contains two tricritical points (A and B) and a triple point 
(C). In the domain between the solid lines BC and CA and the line of circles BCA, FS 
can be overheated. Domains of overcooling of phases FM and N below the solid lines BC 
and CA also exist [8,9]. For example, the domain below the solid line CA and above the 
positive axes is a region of overcooling of the N-phase. 

3. Magnetic susceptibility 

The magnetic susceptibility x — (x/V) = d 2 f(M)/dM 2 per unit volume (V) in a zero 
external magnetic field is obtained in the general form 



where the equilibrium magnetization M should be taken for the respective equilibrium 
phase: (a) M = in the N-phase, (b) M = y/\t\ in FM, and (c) in FS, M is given as the 
maximal nonnegative root of the equation [8,9]. 



We shall compare the known paramagnetic (xp = l/2i;i > 0) and the ferromagnetic 
(Xf = l/4|t|;t < 0) susceptibilities with that corresponding to FS. Note, that both xp 
and xf are divergent on the N-FM phase transition line. While the magnetic susceptibility 
in phases N and FM and the behaviour of the same quantity at the N-FM phase transition 
do not depart from the standard predictions [6] we shall see that the phase transitions 
from FS to N or FM show quite special magnetic properties. As x cannot be analytically 
calculated for the whole domain of stability of FS, we shall consider the close vicinity of 
the N-FS and FM-FS phase transition lines. 

Near the second order phase transition line on the left of the point B (t < — 7 2 /4), the 
magnetization has a smooth behaviour and the magnetic susceptibility does not exhibit 
any singularities (jump or divergence). For t > 7 2 /2, the magnetization is given by 
M = (s- + s+), where 



For r = 0, M = 0, whereas for \^r\ <C (t — 7 2 /2) and r = one may obtain M 
— 7?-/ (2t — 7 2 ) <C 2t. This means that in a close vicinity (r < 0) of r = along the second 
order phase transition line (r = 0,t > 7 2 ) the magnetic susceptibility is well described by 
the paramagnetic law xp = (l/2t). For r < and t — > 7 2 /2, we obtain M = —{^r/2) 1 ^ 
which yields 




(2) 




(3) 






(5) 
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On the phase transition line AC we have [8,9] 

1 



M eq (t) 



2V2 



7 2 - 8t + 7 ( 7 2 + 16£) 



1/2 



1/2 



and, hence, 



7 



16* 

1 + — 

7 2 



1/2' 



(6) 



(7) 



At the tricritical point A this result yields x — 0, whereas at the triple point C with 
coordinates (0, 7 2 /4) we have x(C) = (2/7 2 ). On the line BC we obtain M = 7/2 [8,9] 
and, hence, 



X 



2t + 



r 



(8) 

0. 



At the tricritical point B with coordinates (— 7 2 /4, 7 2 /2) this result yields x 

The comparison of Eqs. (7) and (8) with xp corresponding to the N-phase and with xf 
corresponding to FM shows that the magnetic susceptibility undergoes finite jumps at the 
first order N-FS and FM-FS transitions. The finite jumps Axp and Axf can be easily 
calculated with the help of Axp,f = (x ~ Xp,f) and Eqs. (7) and (8). In particular, the 
jump of x at the triple point C is infinite because of the divergency of xp m the limit 

4. Entropy and specific heat 

The entropy S(T) = (S /V) = -Vd(ffdT) and the specific heat C(T) = (C /V) = 
T(dS/dT) per unit volume V are calculated in a standard way [6]. We are interested 
in the jumps of these quantities on the N-FM, FM-FS, and N-FS transition lines. The 
behaviour of S(T) and C(T) near the N-FM phase transition and near the FM-FS phase 
transition line of second order on the left of the point B (Fig. 1) is known from the 
standard theory of critical phenomena (see, e.g., Ref. [6] and for this reason we focus our 
attention on the phase transitions of type FS-FM and FS-N for (t > — 7 2 /4), i.e., on the 
right of the point B in Fig. 1. 

Using the equations for the order parameters ip and M [8,9] and applying the standard 
procedure for the calculation of S, we obtain the general expression 



S(T) = 



atf 

7J3 



M 



(9) 



The next step is to calculate the entropies Sfs{T) and Sfm of the ordered phases FS and 
FM. Note, that we work with the usual convention Fn = VJ'n = for the free energy of 
the N-phase and, hence, we must set Sn(T) = 0. 

Consider the second order phase transition line (r = 0, t > 7 2 /2). Near this line Sfs(T) is 
a smooth function of T and has no jump but the specific heat Cfs has a jump at T = T s , 
i.e. for r = 0. This jump is given by 



AC FS (T S ) = 



a 2 T, 



1 

2t(T,)/ 7 2 



(10) 



5 



The jump AC FS (T S ) is higher than the usual jump AC(T C ) = T c a 2 /b known from the 
Landau theory of standard second order phase transitions [6]. 

The entropy jump ASac(T) = Sps(T) on the line AC is obtained in the form 



AS AC (T) = -M ( 



eq 



1+ l 1+ 72 



16^ 1/2 ' 



a f 



^M, } , (ID 



where M eq is given by Eq. (6). From (6) and (11), we have AS(t = 7 2 /2) = 0, i.e., AS(T) 
becomes equal to zero at the tricritical point A. Besides we find from (6) and (11) that 
at the triple point C the entropy jump is given by 

On the line BC the entropy jump is defined by AS BC {T) = [S FS (T) - S FM (T)]. We 
obtain 

At the tricritical point B this jump is equal to zero as it should be. The calculation of the 
specific heat jump on the first order phase transition lines AC and BC is redundant for 
two reasons. Firstly, the jump of the specific heat at a first order phase transition differs 
from the entropy by a factor of order of unity. Secondly, in caloric experiments where the 
relevant quantity is the latent heat Q = TAS(T), the specific heat jump can hardly be 
distinguished. 

5. Concluding remarks 

The results obtained in this paper can be used in interpretations of magnetic and caloric 
experiments in various classes of magnetic superconductors with a spin-triplet Cooper 
pairing. Our general consideration is irrespective of the ratio between T s and Tf. But one 
should be aware that depending on the concrete choice of the substance, either T s < Tf 
or T s > Tf may happen. For the new class of ferromagnetic compounds, mentioned in 
Sec. 1, T s <C Tf, and this must be taken in mind when the present theory is compared 
with the experiments performed with these metallic compounds. The condition Tf > T s 
leads to a reduction of the possible phases and phase transitions. For Tf > T s , the 
diagram shown in Fig. 1 is reduced to the domain in the plane (t, r) where the inequality 
r > (a s /ctf)t is fulfilled. Experimental data for the density of states and, hence, for the 
parameters a s and af are needed in order to locate the respective part of the plane (t, r). 
According to the current literature [7] the respective experimental data for the metallic 
compounds enumerated in the Sec. 1 are not yet available. It is, however, obvious that 
the metallic compounds with Tf ^> T s are described by the phase transition lines in the 
fourth quadrant of the (t, r) plane and, therefore, one may observe second or first order 
FM-FS transitions depending on the ratios a s /af and T s /Tf. 

Acknowledgments. Financial support through SCENET (6FP-EC) and Project 10- 
98/2003/JINR (Dubna) are acknowledged. 



6 



REFERENCES 



[1] SAXENA S. S., P. AGARVAL, K. AHILAN, F. M. GROSCHE,, R. K. W. HASEL- 
WIMMER, M. J. STEINER, E. PUGH, I. R. WALKER, S. R. JULIAN, P. MONTHOUX, 
G. G. LONZARICH, A. HUXLEY, I. SHEIKIN, D. BRAITHWAITE, J. FLOUQUET, 
Nature (London) 406, 2000, 587-592. [2] PFLEIDERER C, M. UHLATZ, S. M. HAY- 
DEN, R. VOLLMER, H. V. LOHNEYSEN, N. R. BERHOEFT, C. G. LONZARICH, 
Nature (London) 412, 2001, 58-61. [3] AOKI D., A. HUXLEY, E. RESSOUCHE, D. 
BRAITHWAITE, J. FLOUQUET, J-P. BRISON, E. LHOTEL, C. PAULSEN, Nature 
(London) 413, 2001, 613-616. [4] BLAGOEVA E. J., G. BUSIELLO, L. DE CESARE, 
Y. T. MILLEV, I. RABUFFO, D. I. UZUNOV, Phys. Rev. B42, 1990, 6124-6137. [5] 
UZUNOV D. I. In: Advances in Theoretical Physics, ed. E. Caianiello, Singapore, World 
Scientific, 1990, 96-135. [6] UZUNOV D. L, Theory of Critical Phenomena, Singapore, 
World Scientific, 1993, 1-452. [7] WALKER M. B., K. V. SAMOKHIN, Phys. Rev. Lett. 
88, 2002, 204001-204004. [8] SHOPOVA D. V., D. I. UZUNOV, Phys. Lett. A 313, 
2003, 139-143. [9] SHOPOVA D. V., D. I. UZUNOV, J. Phys. Studies, 7 No 4 (in press). 



7 



